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HEAVING  AND  PITCHING  SHIPMODEL.  WITH  ZERO  FOR.VARD  SPEED  IN  STILL  WATER, 
by  Prof.ir  J.  Gerritt:..a  and  W.  Beukelman. 


Summary. 

Aa  an  extension  of  earlier  work  on  the  distribution  of  the 
hydrodynamic  forces,  acting  on  a  heaving  and  pitching  shipmodel, 
this  report  deals  with  the  case  of  zero  forward  speed. 

Forced  oscillation  teotc  were  carried  out  with  a  segmented 
shipmodel  to  investigate  the  distribution  of  the  forces  along  the 
hull  for  heaving  and  pitching  motions.  The  in-phase  and  quadrature 
components  of  the  vertical  forces  on  each  of  the  seven  sections 
were  measured  as  a  function  of  frequency  and  an  analysis  was  made 
of  their  distribution  along  the  length  of  the  shipmodel. 

The  experimental  results  are  compared  with  Orim'a  values  for 
damping  and  added  mass.  The  comparison  shows  u  satisfactory  agree¬ 
ment  between  theory  and  experiment. 
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1#  Introduction. 


According  to  strip-theory  the  damping  coefficients!  ths 
added  mass  and  the  added  mass-moment  of  inertia  of  a  ship,  per¬ 
forming  harmonic  heaving  and  pitching  motions  in  still  water, 
can  be  calculated  by  integrating  the  corresponding  two-dimen¬ 
sional  quantities  for  each  cross-section  over  the  length  of  the 
ship. 


This  very  simple  procedure  does  not  include  three-dimen¬ 
sional  or  forward  speed  effects  on  the  total  damping  end  added 
mass.  Forced  oscillation  tests  with  a  2. 3  meter  model  of  a  nor¬ 
mal  cargo&hip  have  shown  a  very  small  speed  dependency  of  these 
quantities,  at  least  for  frequencies  of  interest  for  the  usual 
Beekeeping  problems  .  The  speed  range  was  Fn  ■  .15  to  Fn  *  .30. 
The  testa  wer«*  c  ~ried  out  in  the  4.3  meter  wide  towing  tank  of 
the  Delft  Shipbuilding  Laboratory.  Teats  at  zero  forward  speed 
were  not  considered  in  this  tank  because  a  large  wall  influence 
was  expected  from  the  reflection  of  waves,  generated  by  the  forced 
model  motions.  Without  zero  speed  experiments  fora  effects  and 
epeed  effects  could  not  be  separated. 

In  a  later  stage  it  was  decided  to  carry  out  zero  speed  tests 
in  a  wide  basin  to  minimize  wall  interference.  The  results  could 
than  be  used  to  study  form  effect  in  damping  and  added  mass  by 
comparing  the  results  with  calculations  from  strip  theory.  ^The 
ero  speed  experiments  could  also  double  the  range  of  speeds  con- 
cidered). 

This  is  of  interest  for  the  analysis  of  the  damping  cross¬ 
coupling  coefficients  of  the  motion  equations.  These  cross-coup- 
ling  coefficients  depend  quite  strongly  on  the  speed  of  the  ship 
and  this  io  *  *  oortant  for  the  calculation  of  ship  motions  in 
waves.  From  atrip-theory 3  it  follows  that: 


e 


E 


J  N'  x  dx  -  Vm 

L 

I  N*x  dx  ♦  Vm 


(1) 
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whe re : 


•  ia  the  coefficient  of  the  pitch  velocity  tern  in  the  notion 
equation  for  heave, 

E  is  the  coefficient  of  the  heave  velocity  tern  in  the  oquation 
for  pitch, 

N*  ia  the  sectional  damping  coefficient, 

m  is  the  total  added  mass, 

V  ia  the  forward  speed  of  the  ship. 

By  plotting  e  and  E  on  a  base  of  forward  speed  two  straight  lines  cor¬ 
respond  to  each  frequency  of  motion*  The  two  linsa  have  a  common  point 

for  V  =  0,  where:  •«£*  J  N'  xdx  and  their  slope  ia  given  by  ±m* 

L 

The  experiments  have  shown  the  validity  of  equation  (1)  in  the 
epeed  range  Fn  3  0.15*  «30>  The  sero  forward  speed  teats  could  give  ad¬ 
ditional  support  to  the  theory  because  of  the  considerable  extension  of 
the  speed  range* 

Funds  became  available  from  the  Office  of  Naval  Research  by  which 
forced  oscillation  tests  with  a  segmented  model  at  sero  speed  could  be 
carried  out  in  the  Seakeeping  laboratory  of  tha  Netherlands  Ship  Model 
Basin  at  N ageningen. 

In  this  report  the  results  of  the  experiments  are  given.  The  re¬ 
sults  ere  compared  with  calculated  values,  using  Orim'a  two-dimensional 
results  for  damping  and  added  mass  for  ship  like  sections  £2]. 

2.  Experiments. 

The  oscillation  tests  were  carried  out  with  m  2.3  meter  medel  of 
the  Sixty  Series,  having  a  blockcoefficient  Cg»0.7Q.  The  main  dimen¬ 
sions  are  given  in  Table  I.  The  model  is  made  of  polyester,  re-inforced 
with  fibreglass,  and  consists  of  seven  separate  sections  of  equal  length. 
Each  of  the  sections  has  two  end-bulkheads.  The  width  of  the  gap  between 
two  sections  is  one  millimeter.  The  sections  are  not  connected  to  eaoh 
other,  but  they  are  kept  in  position  by  means  of  stiff  strain-gauge  dy¬ 
namometers,  which  are  connected  to  a  longitudinal  steel  box  girder  sbovs 
the  model. 

The  dynamometers  are  sensitive  only  for  forces  psrpendloulsr  to  ths 


bsseline  of  the  model* 


By  means  of  a  Scoth-Yoke  mechanism  a  harmonic  heaving  or  pitch¬ 
ing  motion  can  be  given  to  the  combination  of  the  seven  sections,  which 
form  the  shipmodel.  The  total  forces  on  each  section  could  be  measured 
ae  a  function  of  frequency. 

A  non  segmented  model  of  the  same  form  was  r16o  tested  in  the  same 
conditions  of  frequency  to  compare  the  forces  on  the  whole  model  with 
the  sum  of  the  section  results  *  A  possible  effect  of  the  gaps  between 
the  sections  could  be  detected  in  this  way.  The  arrangement  of  the  tests 
with  the  segmented  model  and  with  the  whole  model  is  given  in  Fig.  1. 

The  mechanical  oscillator  and  the  measuring  system  is  shown  in 
Fig.  2.  In  principle  the  measuring  system  is  similar  to  the  one  des¬ 
cribed  by  Goodman  ^3j:  the  measured  force  signal  is  multiplied  by  costtft 
and  slnU)t  and  after  integration  the  first  harmonics  of  the  in-phase  and 
quadrature  components  can  be  found  without  distortion  due  to  vibration 
noise.  In  some  details  the  electronic  circuit  differs  somewhat  from  the 
description  ir.  CO  .  In  particular  synchro  resolvers  are  U3ed  instead  of 
slne-coslr.e  potentiometer,  because  thev  allow  higher  rotational  speeds* 

The  accuracy  of  tic  instrumentation  proved  to  be  satisfactory  which 
is  important  for  the  determination  of  the  quadrature  components,  which 
are  small  in  comparison  with  the  in-phase  components  of  the  measured 
forces.  Throughout  the  experiments  only  first  harmonica  were  deter¬ 

mined.  It  should  be  noted  that  non-linear  effects  may  be  important  for 
the  sections  at  uhe  bow  and  the  stern  where  the  ship  is  not  wall-sided. 
The  forced  oscillation  tests  were  carried  out  for  frequencies  up  to 
CJ  «  14  rad/sec. 

The  motion  amplitudes  of  the  shipmodel  covered  a  sufficiently  large 
range  to  study  the  linearity  of  the  measured  values.  An  example  of  the 
measured  forces  on  section  2,  when  the  combination  of  the  seven  sections 
performs  a  pitching  motion,  is  given  in  Fig.  3* 

The  position  of  the  model  in  the  tank  is  given  in  Figure  4.  There 
are  beaches  or.  two  adjacent  sides  of  the  tank  and  a  wire  mesh  along  a 
third  side  as  indicated  in  the  figure. 

Each  oscillation  test  was  completed  within  the  time  that  obssrvablt 
reflected  waves  reached  the  model. 
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TABLE  1.  Main  particulars  of  ths  shipmodel. 


Length  between  perpendiculars 

2.2  58 

m 

Length  on  the  waterline 

2.296 

a 

Breadth 

0.322 

a 

Draught 

0,129 

a 

Volume  of  displacement 

0.0657 

§3 

Blockcocfficient 

0.700 

Coefficient  of  mid  length  section 

O.986 

• 

Prismatic  coefficient 

C.710 

..aterplane  area 

0.572 

a2 

Waterplane  coefficient 

0.785 

Longitudinal  moment  of  inertia 

of  waterplane 

0.1685 

v? 

L.C.B.  forward  of  L  /2 

PP 

Centre  of  effort  of  waterplane 

0.011 

m 

after  L  /2 

PP 

0.038 

a 

Froude  number  of  service  speed 

0.20 

3.  Presentation  of  the  results. 

3.1.  //hole  model. 

It  is  assumed  that  the  force  F  and  the  moment  M  acting  on 
a  forced  heaving  and  pitching  shiprodel  in  still  water  can  be 
described  by  the  following  equations: 

Heave : 

(aeP^)a  ♦  b  l  ♦  c  t  »  F  ein(UJt  ♦  d)  , 

"  °  °  }  (2) 

D  E  ♦  E  t  ♦  Ci  *  -M  sinftOt+A)  ' 
o  o  o  s 

Pitch? 

(A*ky2yfV)tUB$eC©  -  M0sinfc*)t«>}f  ) 

d  8  +  e6  +  gtt  -  -  F^in^t  +i ) 

where:  z  is  the  heave  displaceaent  and  0  is  the  pitch  angla. 
o 


6 
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For  a  given  heaving  motion:  z^-z^ainWt,  it  follows  that: 


b  = 


F^ainCt 

z  10 
a 


c  z  -  F  coood 
a  z 

z  CO2 
a 


-  Hz  BlnA 

z  CO 
a 

g  *  ♦  h  co,fi 


-f 


(4) 


s  CO 

a 


Similar  expressions  are  valid  for  the  pitching  motion.  The 
determination  of  the  damping  coefficients  b  and  B  and  the  damping 
cross-coupling  coefficients  e  and  E,  is  straight  forward:  for  a 
given  frequency  these  coefficients  are  proportional  to  the  quadra¬ 
ture  components  of  the  forces  or  moments  for  unit  amplitude  of  mo¬ 
tion.  For  the  determination  of  the  nr  ied  mss,  the  added  nans  mo¬ 
ment  of  inertia,  a  and  A,  and  the  added  naan  cross-coupling  coeffi¬ 
cients  d  and  1)  it  is  necessary  to  know  the  restoring  force  and  mo¬ 
ment  coefficients  c  and  Cf  and  the  statical  cross-coupling  coeffi¬ 
cients  g  and  Q. 

For  zero  speed  of  advance,  these  coefficients  can  be  calculated 
sceily  from  the  known  values  of  the  waterplane  area,  the  longitudinal 
moment  of  inertia  of  the  waterplane  area  and  the  position  of  the  cen¬ 
tre  of  effort  of  the  waterplane  area,  as  given  in  Table  1. 

The  results  for  the  whole  model  are  given  in  Figures  5  and  6. 
These  figures  also  give  the  sum  of  the  results  from  each  section  to 
show  that  the  combined  effect  of  the  slits  between  the  sections  and 
possible  inaccuracies  in  the  measurements  is  very  small. 


3.2.  Results  for  the  sections 


The  cooponents  of  the  forces  on  each  of  the  seven  sections 
were  determined  in  the  sane  way  as  for  the  whole  nodel.  As  only 
the  forces  and  moments  on  the  sections  were  measured  two  equa¬ 
tions  remain  lor  each  section: 

Heave : 

(a*+^V#)z  ♦  b**  ♦c*z  ■  F°  sin(COt+o£.  ) 

1  o  o  o  s 

Pitch: 

(d’^V'x^S^O*  ge«-F*sinfWt  ♦fT*) 

where  (f\7  x ^  is  the  mass  moment  of  the  section  i  with  respect  to 
the  pitching  axis.  The  star  (*)  indicates  the  coefficients  of  the 
sections.  The  section  coefficients  divided  by  the  length  of  the  sec¬ 
tions  give  the  mean  cross-section  coefficients,  thus: 


and  so  on.  Acsuming  that  the  distributions  of  the  cross-sectional 
values  of  the  coefficients:  a',  b'  etcetera,  are  continuous  curves, 
these  distributions  can  be  determined  from  the  seven  mean  cross- 
section  values. 


In  Figure  7  the  distributions  of  the  added  mass  a,  the  damping 
coefficient  b  and  the  cross-coupling  coefficients  d  and  e  are  given 
as  a  function  of  frequency. 

Numerical  values  of  the  section  results,  a*,  b*  etcetera,  are 
summarized  in  Table  2.  In  this  Table  the  results  for  the  whole  model 


are  also  given. 

The  sum  of  the  section  resultw,  as  used  in  Figures  3  and  6  for 
comparison  with  the  whole  nodel  results,  are  determined  according 


to  the  following  equations: 


Zb'  «  b 


Z  d*  =  d 

z  •*  *  • 


X  dx  n  A 
x  dx  ■  B 
x  dx  «  D 
x  dx  »  E 


1 

i 


(6) 


The  values  for  A,  B,  D  and  £  are  given  in  Table  3 
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TABLE  2#  Added  mass,  deaping  coefficient  end  cross-coupling 
eoefficiente  for  the  eectione  and  for  the  whole  aodel. 


U) 

rad/ 

sec. 


sue  of 

whole 

sections 

aodel 

kg  aec2/* 


1.29 

1.20 

1.47 

1.57 

1.68 


kg  aec/m 


2.90 

2,78 

0,68 

1.75 

1,17 

0.21 

1.26 

0.88 

0.18 

1.22 

0.85 

0.11 

1.36 

0.85 

0.12 

1.63 

0.90 

0.14 

1.73 

1.04 

0.17 

2 

0 

0 

• 

r 

1.40 

3.95 

3.90 

3.20 

4 

2.08 

4.48 

5.18 

5.93 

5.75 

6 

1.98 

4,23 

5.10 

5.15 

4.62 

8 

2.15 

3.80 

2.93 

2,13 

2.81 

10 

2.18 

3.06 

1.72 

1.28 

1.31 

12 

1.98 

2,22 

1.32 

0.97 

0.64 

14 

1.54 

1.78 

1.51 

1.03 

0.69 

- 

kg  sec" 


.06 

.03 

0.01 

0.00 

0.00 

0.00 

0.00 


-0.01 

-0.01 

-0.01 

♦0.09 


o.8i 

0.41 

0.30 

0.30 


•  33 
i.3o 
0.56 
0.45 
0.45 

0.51 

0.53 


0.74 

0.56 

0.42 

0.30 


5.51 

6.02 

6.58 

7.14 


17.65 

29.72 

25.80 

17.09 

11.89 

8.96 

8.33 


1.19 

Q 

2.10 

■  1 

1.62 

0.72 

0.95 

0.54 

O.69 

0.45 

-0.63 

O.38 

—0.46 

-2.09 

-2.53 

-2.57 

-2.21 

-1.94 


18,80 

8.61 

5.84 

5.43 

5.86 

6.48 

6.82 


10.25 

29.17 

25.03 

18.00 

12.31 

9.43 

7.87 


♦0.13 
-O.69 
-1.81 
-2.46 
-2.56 
-2.08 
.95 


TABLE  3.  Added  mas*  moment  of  inertia,  damping  coeeficient  and 
cross-coupling  coefficient-*  xor  the  whole  model. 


CO 

2 

A  kgcsec 

B  kg  a  sec 

2 

D  kg  sec 

L  kg  sec 

rad/ 

sum  of 

3  UK  Of 

whole 

sum  of 

whole 

sum  of 

whole 

sec. 

sections 

sections 

mof’.cl 

sections 

model 

sections 

model 

2 

3.77 

0.44 

• 

-1.28 

• 

♦0.55 

4 

2.26 

2.34 

3.72 

3.72 

-oe58 

-0,41 

-0,52 

-0.63 

6 

1-54 

1.49 

6.94 

6.90 

-0.19 

-0.18 

-1.63 

-1.83 

8 

1.18 

1.13 

5.83 

5.65 

-0.08 

-0.06 

-2.42 

-2.49 

m  i 

1.13 

1.09 

4.50 

4.47 

♦0.06 

♦0.02 

-2.58 

-2.64 

1.18 

msm 

3.53 

3.29 

fO.08 

♦0.03 

-2.46 

-2.32 

BB 

- 

KB 

- 

2.84 

• 

♦0.04 

am 

-2.07 

For  A  and  B  the  6un  of  the  section  results  were  obtained  by 
using  equation  (6). 


4 .  Analysis  of  the  results. 

The  experimental  values  of  the  coefficients  of  the  motion  equa^> 
tions  will  now  be  analysed  by  using  the  strip-theory.  For  a  descrip¬ 
tion  of  the  strip-theory  reference  is  made  to  £lj  •  In  the  special  case 
of  zero  forward  speed  the  coefficients  of  the  motion  equations,  accor¬ 
ding  to  the  strip-theory,  have  a  very  simple  fora,  as  shown  in  Table  4, 
where: 


a 

A 


.  .  g  .  Coefficients  of  the  motion  equations 
.  .  GJ  (hydromechanical  part). 

a*  Added  mass  of  a  cross-section  (zero  speed) 

N'  Damping  coefficient  of  a  cross-section  (zero  speed). 

Aw  .Vaterplane  area. 

Sw  Statical  moment  of  waterplane. 

Iw  Longitudinal  moment  of  inertia  of  waterplane. 
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TABLE  4.  Coefficients  of  the  motion  equations  according  to 
strip-theory  for  zero  forward  speed. 


Table  4  shows  that  each  of  the  coefficients  cnn  be  calculated 
by  integrating  the  corresponding  croas-sectional  values  over  tho 
length  of  the  ship. 

This  procedure  has  been  carried  out  for  c he  c.odol  under  consi¬ 
deration.  Grim's  cross-sectional  values  for  the  damping  coefficients 
N'  and  the  added  mass  m'  were  used. 

In  Figure  7  the  calculated  distributions  of  the  added  mass  and 
the  damping  coefficients  are  shown  in  comparison  with  the  experimen¬ 
tal  values.  In  genoral  there  is  satisfactory  agreement  between  theo¬ 
ry  and  experiment.  The  calculated  added  mass  at  the  lowest  frequen¬ 
cies  is  smaller  than  indicated  by  the  experiments,  and  some  of  the 
distributions  show  a  slight  shift  in  longitudinal  direction  when 
compared  with  the  measured  values. 

Also  the  values  of  the  total  added  mass,  the  damning  coeffi¬ 
cients  und  the  cross-coupling  coefficients  are  also  compared  in  the 
Figures  8  and  9» 

The  numerical  values  are  summarized  in  Table  5. 
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TABLE  5%  Conparia on  of  calculated  and  measured  coefficients 


-  12  - 


la  the  region  of  resonance  (for  heave  CO  •  6.9  rad/eec, 

for  pitch  CJ  =  7.0  rad/sec)  the  agreement  between  theory  and 
n 

experiment  is  excellent. 

For  frequencies  lowwr  than  5  rad/sec  the  calculated  pitch 
damping  coeff. cient  and  the  added  mass  moment  of  inertia  show  a 
rather  large  difference  between  calculation  and  experiment.  The 
heave  and  pitch  damping  coefficients  for  10  rad/sec  are 

smaller  than  the  measured  values. 

Although  only  first  harmonic  component c  of  the  hydromechani¬ 
cal  forces  were  measured,  there  is  no  garantee  that  the  90  degrees 
out  of  phase  components  correspond  exclusively  to  wave  damping. 
Particularly  at  high  frequencies  damping  of  viscous  origin  could 
be  present  and  the  first  harmonic  content  of  such  p  non  wave  making 
damping  could  influence  the  final  result.  Because  the  calculation 
only  considers  pure  wave  damping  this  could  be  a  reason  for  the  ob¬ 
served  differences. 

For  practical  purposes,  ouch  as  the  calculation  of  ship  mo¬ 
tions  in  waves,  tho  differences  as  shown  in  the  Figures  8  and  9  are 
not  of  much  interest. 

As  3  final  comparison  between  experiment  and  atrip-theory,  the 
calculated  and  measured  damping  cross-coupling  coefficients  e  and  E 
at  zero  forward  speed,  are  plotted  in  Figure  10,  which  is  taken  from 
reference^! J.  In  Figure  10  it  is  shown  that  e  and  E  vary  linearly 
with  forward  speed  at  constant  frequency  of  motion.  The  range  of  for¬ 
ward  speed  were  this  is  valid  now  includes  zero  speed  and  the  calcu¬ 
lated  damping  croac-coupling  coefficients  are  in  agreement  with  the 
measurements. 
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7.  List  of  symbols. 


a  . 
A  . 

a*. 

A*. 

e'. 

A’. 


gi  Coefficients  of  the  motion  equations. 
qJ  (hydromechanical  part). 

*  \The  same  for  a  section  of  the  ship. 
0*1 


)The  same  for  a  cross-section  of  the  ship. 


Aw  Waterplane  area. 

Cg  Blockcoef ficient. 

Fn  Froude  number. 

F£,  Fg  Amplitude  of  vertical  force  on  a  heaving  or  pitching 
ship. 

g  Acceleration  of  gravity. 
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Iv  Longitudinal  moment  of  inertia  of  waterplane 

Longitudinal  radius  of  inertia  of  the  ship. 

L^  Length  between  perpendiculars. 

Mg |  Mg  Amplitude  of  moment  on  a  heaving  or  pitching  ship, 
m'  Added  mass  of  a  cross-section  (aero  speed), 

m  Added  mess  of  a  ship  (heave). 

N'  Damping  coefficient  of  a  croas-aection  (zero  speed). 
Sw  Statical  moment  of  waterplane. 
t  Time. 

V  Forward  speed  of  ship. 

x«  y,  z  Right  hand  coordinate  system,  fixed  to  the  ship. 


x  ,y  ,  z 
owo  o 


0 

9. 

p 

O) 

V 

V 


Right  hand  coordinate  system,  fixed  in  space. 

Vertical  desplacement  of  ship. 

Heave  amplitude. 

Horizontal  distance  of  centre  of  gravity  of  a  section  to  the 
pitching  axis. 

Phase  angles. 

Pitch  angle. 

Pitch  amplitude. 

Density  of  water. 

Circular  frequency. 

Volume  of  displacement  of  ship. 

Volume  of  displacement  of  section. 
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ARRANGEMENT  OF  OSCILLATION  TESTS 


PITCH  -60 

SECTION  2  (AFT) 


FIG.3  COMPONENTS  OF  FORCE  ON  SECTION  2.  PITCHING  MOTION 


INSTRUMENTS 


FIG.  4  POSITION  OF  MODEL  IN  SEAKEEPING  BASIN 
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FIG.  5  EX 


d-— COUPLING  COEFFICIENT  bgs«c  A — -ADDED  MASS  MOMENT  OF  INERTIA  kgmicc 
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FIG  8  COMPARISON  OF  THE  CALCULATED  VALUES  OFo.b.D  AND  E 
WITH  THE  EXPERIMENTS 


COMPARISON  OF  THE  CALCULATED  VALUES  OF  A,  B  ,d  AND  • 
WITH  THE  EXPERIMENTS 


